Quantum signatures of breather-breather interactions 
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The spectrum of the Quantum Discrete Nonlinear Schrodinger equation on a periodic ID lattice 
shows some interesting detailed band structure which may be interpreted as the quantum signature 
of a two-breather interaction in the classical case. We show that this fine structure can be interpreted 
using degenerate perturbation theory. 

PACS numbers: 63.20.Pw 



I. INTRODUCTION 

The localization and transport of energy in lattices 
by intrinsic localized modes or discrete breathers, has 
recently been the subject of intense theoretical and ex- 
perimental investigation ( see and references therein). 
Corresponding results on the quantum equivalent of these 
states are less numerous, ci. Hi for some theoretical 
results and Q for some experimental work. Studies of 
quantum modes on small lattices are of increasing inter- 
est for quantum devices based on quantum dots (c.f. 0) 
and for studies of photonic crystals. 

We present some results on breather bands in small one 
dimensional lattices with a small number of quanta. We 
study a periodic lattice with / sites containing bosons, 
described by the quantum version of the discrete nonlin- 
ear Schrodinger equation (DNLS) which will be denoted 
by QDNLS. It has the Hamiltonian 
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which conserves the number of quanta n. This model, 
also known as the Bose-Hubbard model, is presently used 
to investigate cold bosonic atoms in optical lattices Q . 

Assuming that the anharmonic parameter 71 > is 
stronger than the intersite coupling, the eigenvalues of 
l|T)l. plotted as a function of the wave number k, sepa- 
rate out into distinct bands 0- The eigenstates of each 
band are dominated by states in which the bosons have 
clumped together, two or more on one site. For example, 
with n = 4, the lowest band is, to a good approximation, 
a linear combination of states with 4 bosons on site i and 
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FIG. 1: Detail of the eigenvalue spectrum for the QDNLS 
periodic lattice Q, n = 4, / = 19, 71 — 10, e — 0.5. 



no bosons elsewhere. The next lowest band is mostly 
composed of states with 3 quanta on one site and another 
quanta elsewhere. The third band is mostly composed of 
states with 2 quanta on one site and 2 quanta on a sep- 
arate site. We will refer to this band as a {2, 2} band in 
an obvious notation. This band is of great interest since 
it represents the simplest case of a band describing two 
"composite" particles interacting with each other. Our 
letter is devoted to the fine structure of this and similar 
bands such as the {4, 2} and {3, 3} bands in the n = 6 
case. 

The fine structure of the {2, 2} band, see Fig.^ shows 
the eigenvalues (crosses) in the n = 4 case. The solid lines 
show the results of the theoretical calculations (described 
below) in the asymptotic limit /—►(». We stress again 
that this picture shows only one of the bands in the n = 4 
spectrum. Its fine detail is revealed as a "continuum" 
band (in the / — > 00 limit), plus a single fc-dependent 
"line band" . Examination of the corresponding eigen- 



2 



vectors shows that the "line band" is mostly composed 
of states where the two sites each with two quanta lie 
occupy adjacent sites, whereas in the "continuum" band 
these sites are separated by one or more vacant sites. 



II. GENERALIZED HAMILTONIAN 

The standard QDNLS Hamiltonian is thermodynami- 
cally unstable, i.e. the energy of the ground state goes to 
— oo as the number of particles goes to +00. Moreover, 
some bands can be degenerate with others, for example 
when n = 6, the {3, 3} band overlaps with the {4, 1, 1} 
band. In some cases it would be convenient for the bands 
under study to be the lowest energy bands of the system. 

For these reasons we consider also the following gener- 
alized QDNLS Hamiltonian 



H 2 — Hi 
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The 72 > term is a saturation term which discourages 
too many bosons occupying the same site. A similar term 
appears also in nonlinear optics [3 and in cold bosonic 
atoms trapped in optical lattices |8( where it describes 
three body interactions. By varying the values of 71 and 
72 it is possible to "shuffle" the order of bands in the 
(E, k) diagram. 

At zero coupling (e = 0), we denote by E s the respec- 
tive energies of s bosons on the same site, E\ =0, E2 = 
-271, E 3 = -6(71 - 72), E 4 = -12(71 - 272), etc. 
Thus in the n = 4 case the {4} band has energy E4, the 
{31} band has energy E 3 , the {22} band energy 2E 2l etc. 
Providing 71 < 372, the bottom of the {22} band is the 
ground state of the system. 



III. PERTURBATION THEORY 

To describe the components of the quantum states we 
use a position state representation, where for example the 
state 1 020020 . . .) represents a state with two quanta at 
the lattice position 2, two quanta at the lattice position 
5, and zero quanta elsewhere. In view of the periodic na- 
ture of the lattice, we can generate an equivalence class 
of states by applying the translation operator (with pe- 
riodic boundary condition) to one of these states. Wc 
will refer to these classes by ordering them such that the 
leftmost number is the largest, so for example |3000 . . .) 
is shorthand for |3000 . . .), |0300 . . .), |0030 . . .), etc. For 
further conciseness we will truncate all trailing zeros, so 
the above class becomes |3). The set of all the classes 
containing |22), |202), 1 2002) , etc. is referred to as the 
{2,2} band. 

At zero coupling, all the {2,2} states |22), 
1 202), 1 2002), . . . are degenerate. We use degenerate per- 
turbation theory to obtain both eigenvalues and eigen- 
states for e ^ 0. For the sake of simplicity, we consider 



an odd number of sites / = 2<r+ 1. Bloch waves of {2, 2} 
states can be written (in the notation of 0) 



where 



no = E c ^ 
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Here T is the translation operator, r = e lfc , and the crys- 
tal momentum k = 2nl/ f where I £ {— cr, . . . , er}. 

Using © and standard Brillouin-Wigner perturbation 
theory up to second order in e we obtain 



(E-2E 2 ) Cj =J2(^\V\ipj')cf + 

(^\V\^\V\^ r ) 



(5) 
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2E 2 -E 



where V is the hopping term in the Hamiltonian, and 
is any state not in the {2, 2} subspace. E is the energy 
corresponding to in the uncoupled limit (e = 0). 

It is obvious that the first sum of is zero as V does 
not link any of the |^ s ) to each other. In evaluating the 
second sum, the (•|V|*) terms are mostly zero except for 
the following 

(22|V]112) = e 4 V 2 e = e 2fe (22|^|211), 
(22|U|13) = V6e = (22|V|31), 
(202|V|112) = V2e = (2Q2|V|211), 
(202|V|1102) = e*V2e = e ,;fc (202|U|2011), . . . 
(22) 

We now define the matrix H - ■/ by 
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Then 
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where I m is the m x m unity matrix and 

3(^2-2^1) 1 3 72 - 4 7 i 



2E 2 — E 3 — Ei 
k 



1 



7i - 372 



k = e l 2 cos ( — ] and p = cos ok . 
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FIG. 2: Detail of the eigenvalue spectrum for the generalized 
QDNLS periodic lattice ©, n = 4, / = 19, 71 = 10,72 = 
7.5, e = 0.5. 
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FIG. 3: Detail of the eigenvalue spectrum for the generalized 
QDNLS periodic lattice ©, n = 6, / = 11, 71 = 30,72 = 
0,e = 0.5. 



The structure of the matrix J7J) is very similar to the 
two-quanta case described in The first term repre- 
sents a global shift of the {2, 2}-band, whereas the "im- 
purity" r will in general be responsible for a splitting of 
the states |22) from the rest of the band (|202), |2002), 
etc.). In this respect, the |22)'s can be seen as "bound 
states of doublets" within the {2, 2}-band. It is separated 
from the continuum band because in addition to linking 
with states such as 1 2 1 1) , it also links to states such as 
|31). This explains the {2, 2}-band fine structure. 

In case the number of sites tends to infinity, Q may 
be diagonalised exactly to yield 



120.04- 



120.00- 



119.96 




fA(2 + r + ^F) iff |r| > cos (I) 

E = 2E 2 +l 

{ E s f 2E (1 + cos (|) cos 9) where e (0, tt) 

(10) 

Depending on the various values of 71,72, and e, we 
can get the line band completely above or below the con- 
tinuum band, or partially merged with the continuum 
band. Fig [5] shows a case where the line part is below 
the continuum part of the {2, 2} band and represents the 
ground state of the n = 4 sector. 

To provide some insight in the way the two groups of 
two quanta interact within a 1 22) state, we compare its ef- 
fective mass ?7i22 with twice the effective mass m\ of a sin- 
gle |2) state. From 2] and ifTnjl . we obtain m\ ~ 7i/(2e 2 ) 
as e — > and m^/f^) = T (see ©). Depending on 
the ratio 72/71, T and thus are either positive or 
negative. The same phenomenon is observed for bright 
solitons in systems of Bose-Einstein condensates in opti- 
cal lattices for which JTJ may be seen as a tight-binding 
limit (see [Tol | and references therein). 



FIG. 4: Detail of the eigenvalue spectrum for lj2"{. n — 6. 

/ = 11, 71 = 10, 72 = 20, e = 0.5. 



IV. THE n = 6 CASE 

As a further example we consider the n = 6 case. Now 
there are two bands describing an interaction of two an- 
harmonic states: the {4,2} and the {3,3} bands. Fig [31 
shows a n = 6, {4, 2} example, again the crosses repre- 
sent numerically exact solutions in the / = 11 case, and 
the lines represent perturbation theory calculations. In 
this case we have a "continuum band" which shows very 
weak fc-dependence, plus two "line bands" , one above and 
one below. With other choices of 71 and 72 we can move 
one or both of the "line bands" into the continuum. 

Fig^Jshows the corresponding {3, 3} case. In this case 
there appears to be only two "line bands" , but a closer 
examination reveals the upper "line" is in fact a contin- 
uum band with the degeneracy split at an 0(5/7) 3 level. 
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A. Perturbation theory, n = 6, {4, 2} band 

Proceeding as in the {2, 2} case, we now obtain 

P\ 



ff< 4a >=De a J 2ff - — 
7i 



(T 1 
' 1 1 



1 1 
1 Tj 



(11) 



r responsible for a splitting of the states 1 33) from the 
rest of the band. In this respect, the |33)'s can be seen as 
"bound states of doublets" within the {3, 3}-band. Note 
again that none of the elements of the matrix above con- 
tains the wave vector k. At this order of perturbation 
theory, the bands are then flat. Moreover, the degen- 
eracy of the |3 • • • 3)'s (i.e. two 3-quanta breathers sepa- 
rated by one or more empty sites) has still not been lifted 
although it would to next order. 



where 



D _ 2 571 ~ 972 . r 2 



4 7l 2 - 27 7 | 



P : 



3 7i(7i - 372) 
67ie lfc 



3 (71 - 372) (71 - 672) 



71 - 672 



An analysis of the eigenvalues of this matrix in the / — > 
00 limit is somewhat messy but explicit results can be 
found: details will be published elsewhere. Essentially 
the results depend on whether two rational functions of 
71/72 are less than or greater than 1 in modulus. This 
gives the two "line bands" in Fig EI The corresponding 
eigenstates are essentially symmetric and antisymmetric 
combinations of the Bloch waves made from |42) and 1 24) . 
The "continuum band" is given by the formula 



E c = 24 72 - 14 7 i 



2e2 

7i 



571 - 972 



972 — 371 



cos 8 



(12) 

where 9 € (0,7r). Note that, up to order 2 in e, the ener- 
gies of this band do not depend on the crystal momentum 
k. in accord with the numerical results. 



B. Perturbation theory, n = 6, {3, 3} band 

In this case the correction to the zeroth order states is 
given by 



#(33) 



Gc 2 



372 - 271 



( U + M) 



(13) 



where Mi,i = T = 5(272 - 71)/ (71 - 672) and M tJ = 
otherwise. The diagonal matrix (|13|l has an "impurity" 



V. DISCUSSION 

It is possible to generalize the above results to discuss 
{to, £} bands when m + I = n > 6. The case m = £ is 
similar to the {3, 3} case and the case to > i > 1 follows 
the {4, 2} case discussed above. Details will be given 
elsewhere. 

The results known previously for bound states repre- 
senting a group of n bosons located on the same site 
have been extended here to higher order states represent- 
ing two interacting groups of bosons. This opens up the 
possibity of a study of the collision process between these 
composite particles, that is a quantum breather collision. 
A similar study has already been done for the continu- 
ous version of ifTfl. i.e. the integrable quantum nonlinear 
Schrodinger equation (QNSE), well known in nonlinear 
optics |TH and quantum field theory |T^ . However, the 
bound states of groups of quanta which we have described 
here within the {to, £} bands do not exist in QNSE. Their 
appearence in these nonintegrable discrete systems is thus 
expected to affect the collision process between two quan- 
tum breathers and to bring new features in comparison to 
the quantum soliton collisions described in [Ij. Results 
of this investigation will be reported elsewhere. 
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